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Color Transfer in Associated Heavy-Quarkonium Production 
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We study the production of heavy quarkonium in association with an additional heavy pair. We 
argue that important contributions may come from phase space regions where three heavy fermions 
are separated by relative velocities much lower than the speed of light, and to which standard 
effective field theories do not apply. In this region, infrared sensitive color exchange is specific to the 
presence of the unpaired (anti) quark. This effect vanishes as the motion of the additional particle 
becomes relativistic with respect to the pair, and is completely absent for massless quarks and gluons 
in the final state. 



In heavy-quarkonium production, the formation of 
the heavy-quark pair can be perturbative. The sub- 
sequent evolution of the quark pair is conventionally 
treated in the language of effective theories, most no- 
tably non-relativistic QCD (NRQCD) Q. In NRQCD, 
non-perturbative dynamics is organized through a factor- 
ization characterized by a joint expansion in the strong 
coupling, a s and v, the relative velocity. This expansion 
is well-justified for heavy-quarkonium decays, although it 
still lacks a compelling proof for production processes. It 
is clear, however, that NRQCD cannot apply directly to 
reactions in which three heavy particles, two quarks and 
an antiquark for example, are produced close together in 
phase space, simply because there is no unique choice of 
relative velocity v. 

Although the application of NRQCD to production 
processes has had many successes [H there remain 
well-known data anomalies, particularly in connection 
with the polarization observations at the Tevatron Q 
and the size of the associated production cross section 
for J/ip with ccpairs, as seen by the BELLE and Babar 
collaborations [5]. Regarding the latter, the rate for J/ip 
associated production is unexpectedly large compared to 
J/tp with light hadrons. The importance of associated 
production within the formalism of NRQCD has recently 
been explored at leading order (LO) for hadroproduc- 
tion in Ref. @ and at next-to-leading order (NLO) for 
e + e~ annihilation in Ref. Q. We will see, however, that 
NRQCD may not be the complete story for these pro- 
cesses. 




FIG. 1: Representative lowest order diagrams for a virtual 
photon decaying into two heavy quark pairs, for fixed quark 
momentum Pi and antiquark momentum P%. The pair Pi 
and Pi will be associated with a bound state. 



Diagrams for e + e _ — > QQ + QQ are shown in Fig. [T] 
Because we will study final states with two pairs, we will 
find it useful to refer to the pair that forms the bound 
state as the "active" pair. To be specific, we will denote 
the quark and antiquark momenta of the active pair as 
Pi and P2, respectively, with total and relative momenta 
and rest-frame relative velocity given by 
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Associated pair production for J/tf) and related 
quarkonium states provides access to a novel kinematic 
region, where the heavy quarkonium is produced close 
in phase space to open heavy flavor. The presence of 
the additional quark, itself a source of color, could influ- 
ence hadronization. In addition, studying this configura- 
tion sheds light on NRQCD factorization H for processes 
without associated heavy quark pairs, and we will make 
contact with our previous work on that important ques- 
tion We begin with a discussion of matching and 
infrared poles in dimensional regularization for NLO cor- 
rections to associated production of heavy quarks Q, of 
mass to, in e + e~ annihilation. 



where 2E* is the total center of mass energy of the active 
pair, and where in the second expression for v 2 we give 
the invariant form. In the diagrams of Fig. [TJ we refer to 
the additional heavy particle that is closest in phase space 
to the active pair as the "spectator" quark (or antiquark) . 

The LO NRQCD cross sections for the amplitudes of 
Fig. [T] are given by 

da e+e -^ H+x {pH) = d K+e-^QQ[n]+x(PH) P") , 
n 

(2) 

where da e + e -_,QQ^ +x {pH ) is an NRQCD coefficient 
function, projected onto heavy quark pair state QQ[n]. 
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The corresponding NRQCD operators are of the form [l| 
0^(0) = E X t (0)«„^(0) \N,H) (N,H\ ^t( )«; x (Q), 

AT 

(3) 

with K n and n' n projections for color (octet or singlet) 
and spin, and where we suppress gauge links [9(. At 
LO, NRQCD singlet configurations are present only in 
diagrams where the active quark and antiquark arise from 
different lines, as in the first graph on the right-hand side 
of Fig. [T] In the second graph, the active pair is purely 
octet, and is expected to contribute at a much lower level 
to the cross section. 

For NRQCD factorization to be predictive, it is neces- 
sary that the coefficient function, da e + e -_ >QQ[ n ]+x(pH) 
be infrared safe at higher orders. Not only infrared sin- 
gularities, but also singular dependence in any of the rel- 
ative velocities must be absorbed into the long-distance 
matrix elements of the effective theory, a process referred 
to as matching. At LO, all gluons are off-shell by at 
least a multiple of m, and the full LO amplitudes may 
be absorbed into coefficient functions. Nontrivial issues 
of factorization and matching first arise at NLO. 

It is a major undertaking to compute the full NLO cor- 
rections to the singlet NRQCD cross section [7|, and all 
the more so for the full four-particle cross section at finite 
relative velocities. Nevertheless, it is relatively straight- 
forward to check the self-consistency of the NRQCD fac- 
torization at NLO by verifying the cancellation of in- 
frared poles. 

As noted above, when more than two relative velocities 
are small, the standard NRQCD power counting does 
not apply directly, simply because there are now three 
dimensionless scales, Pij = ^/l — 4m 2 / (Pj + Pj) 2 , i ^ 
j = 1,2,3, instead of the single v = (3u defined by the 
active pair. Thus, where there was once one "soft" region 
mv and one "ultrasoft" region mv 2 , there are now three 
of each. Let us see how these considerations emerge in 
the NLO diagrams of Fig. O 

We stress that the region of small relative velocity for 
three of the four particles in double heavy-pair produc- 
tion is leading-power in the overall center-of-mass energy, 
y/s. For example, the two diagrams on the right-hand 
side of Fig. Q] contribute to the fragmentation of a heavy 
quark into the quark plus a pair. Such contributions also 
occur as part of the high-pT cross section for associated 
production in hadronic collisions, as studied in Ref. 



It is worth noting as well in this context that in such 
fragmentation- like contributions, the natural scale for a s 
is m, while for the corresponding LO contributions in 
which a light pair recoils from a pair of heavy quarks, the 
natural scale for a s is ^/s. This is an additional source 
of relative enhancement for associated production, com- 
pared to J/ip cross sections with light quanta only. 

As discussed extensively in [9j, factorization requires 
the cancellation of all infrared gluons that are not "topo- 
logically factorized" . Only such contributions will match 
the single-pair matrix elements of NRQCD. We may 
think of a topologically factorized singularity as one that 
can either be absorbed into a quarkonium matrix ele- 
ment or that does not involve either of the active heavy 
quarks. An example of the former is the first diagram in 
Fig. [U and of the latter, the second. After a sum over 
real gluon emission at fixed color for the active quark 
pair, exchanges between spectators cancel. 

The third and fourth diagrams of Fig. [51 however, are 
not topologically factorized, and include the transfer of 
color between the active pair produced at short distances 
and the spectator. A pair produced as an octet in the 
hard scattering may end up as a singlet in the final state. 
Let us see how this mechanism enters the NLO correc- 
tions, and why it is not a problem except in the region 
where the spectator comes close to the active pair. 
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FIG. 2: Diagrams for one-loop virtual infrared corrections, 
where the blob represents a complete set of lowest order Feyn- 
man diagrams, as shown in Fig. [TJ 



When fiij <C 1, the infrared pole behavior of any of the 
exchanges shown in Fig. [2] is given up to a sign by 
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with corrections that are finite as £ = 2— Z?/2 ^ or that vanish as higher powers of /3y . We will concentrate on the 
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imaginary term, which is singular for both vanishing (3 
and e. Consider now the imaginary contributions of Eq. 
(J4j) to the NLO ampliude A, traced over the colors of the 
active pair to enforce a singlet configuration in the final 
state. Denoting by Aij the amplitude with an exchange 
between Pi and Pj, we have, suppressing common color 
factors in the quark and antiquark couplings, 



tator quark (of momentum P3) and P/2 [13 
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where A^ represents the LO amplitude. From this soft 
gluon correction, the third and fourth diagrams in Fig. 
[2] inherit color singlet configurations for the pair Pi , P2 
in the final state at NLO. We should mention that this 
contribution is subleading in the number of colors iV c . 
Indeed, it constitutes an alternative mechanism to the 
"preconfinement" scenario, identified long ago in Ref. 

and built into most hadronization models in event 
generators 

The essential feature of Eq ([5]) is the presence of poles 
in the relative velocities P13 and /?23- These poles cancel 
identically for Pi ■ P 3 = P 2 ■ P 3 or for P 3 2 = 0. In the 
former case, the relative velocity of the spectator quark 
and the active quark equals the relative velocity of the 
spectator and the active antiquark. Here, soft gluons 
emitted by the spectator cannot resolve the charges of 
a perfectly co-moving pair in a singlet color state. Cor- 
respondingly, when the spectator is lightlike, its relative 
velocity to either active line is unity, and the two terms 
cancel. Thus for pair production in the presence of mass- 
less quanta only (light quark pairs and gluons) the color 
transfer mechanism is completely absent. This is an es- 
sential ingredient in the verification of NRQCD factor- 
ization at NNLO in Refs. Evidently, a light-like 
spectator can neither resolve nor produce a color singlet 
even when the quark and antiquark are not co-moving. 
We see here, however, that such a cancellation requires 
that any heavy spectator be moving relativistically with 
respect to the active pair. For nonrelativistic motion, 
there is a nontrivial color transfer from the triplet spec- 
tator quark, that catalyzes the transition between octet 
and singlet. 

If, as in Ref. [7], the active quark relative velocity v 
appears only as a regulator for active-active gluon ex- 
change, nonfactorizing singularities never occur. At the 
same time, at NLO, this cancellation is not exact in the 
low-velocity region for the spectator. The color trans- 
fer contribution is neither topologically factorized nor in- 
frared finite, and cannot be absorbed into the coefficient 
function for the singlet production amplitude. Although 
this effect is purely imaginary at NLO, it will have a 
physical role at higher orders. 

Since the infrared correction to the amplitude, A of ([5]) 
vanishes at vanishing active-pair relative momentum, it is 
natural to expand A in v at fixed relative pair-spectator 
velocity, which we will refer to as (3s, between the spec- 
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with Pq = (2m, 0) and Qs'Pq = 0. Expanding in the rela- 
tive momentum of the active pair in the spirit of NRQCD, 
we find 
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where in the second relation 0s is the angle between the 
vectors qs and q. The presence of the cosine, which is odd 
under the interchange of the members of the active pair, 
is an expression of the dipole nature of this contribution 
to the amplitude. 

Although the expansion in Eq. |(7J) is applicable in only 
a limited region of phase space, it is a region in which the 
active pair and spectator quark form a system of invari- 
ant mass of order (3m) 2 . In this kinematic region, the 
squares of NLO corrections to LO diagrams like those on 
the right-hand side of Fig. [1] contribute to the cross sec- 
tion at leading power in the overall center-of-mass energy. 
At NNLO, the squares of the imaginary poles of Eq. ((4]) 
produce a^(l/e 2 ) corrections to the cross section. All 
other diagrams contribute only at the single-pole level 
[HI]. We note again, that in accordance with the results 
of Ref. Q, this contribution vanishes in the limit m — > 0, 
and so is consistent with NRQCD factorization at large 

Pt- 
Iti summary, we have identified a new, leading-power 
"color-transfer" infrared-divergent contribution to the 
formation of color-singlet pairs, which is specific to as- 
sociated production. This contribution must be added 
to that of the familiar color-octet mechanism .1]. The 
infrared divergence signals that color transfer cannot be 
absorbed into a coefficient function, while its dependence 
on (3s signals that it cannot be absorbed into matrix el- 
ements in an effective field theory based on the dynam- 
ics of a single pair. This mechanism is non-leading in 
the number of colors, although that by itself should not 
eliminate it from consideration. 

It is natural to try and estimate the magnitude of this 
new effect. A plausible estimate can be given by elabo- 
rating on the analogy just observed between the normal 
color octet mechanism and color transfer. In both cases, 
an S-wave pair is transformed to P-wave by the gluonic 
field strength contracted with the pair's color dipole. In 
the octet mechanism, the field couples to the polarization 
of an emitted gluon, while in color transfer it is coupled 
to the Coulomb field of the spectator. To return to an S- 
wave state, we will assume a /3s-independent v 2 suppres- 
sion, following NRQCD power counting [l|, |2j. At fixed 
(3s, let us assume that the nonperturbative formation of 
quarkonium state, H is similar for these two transitions. 
We then square the /3s-dependent kinematic factor of Eq. 
([7]) , and replace the v -dependence by the NRQCD matrix 
element ( 3 S$). In the total cross section at fixed (3s, such 
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a quantity would be additive to the color singlet contri- 
bution, which dominates in the NRQCD treatment. We 



that is, the standard color singlet, plus a color-octet con- 
tribution enhanced by the factor with (3 s the ve- 
locity of the closest heavy spectator, Q 1 = Q,Q. By 
numerical evaluation, we have verified [l3| that the short- 
distance cross section, da e + e ^ ^qq^s 8 ]+qi [p s ){ph) is also 
larger than the corresponding singlet contribution. At 
f3s ~ 0.3, our estimate of color transfer would then be- 
come competitive with the color singlet cross section at 
the same (3s- Of course, for (3s < v an expansion in v/ (3s 
loses meaning, and we do not offer Eq. © as a realistic 
estimate, especially for J/ip, with v ~ 0.5. We do present 
it, however, as an indication of the possible importance 
of the new mechanism. 

A signature of the color-transfer process would clearly 
be an enhancement in any open heavy flavor at low rela- 
tive velocity to quarkonia, in hadronic scattering as well 
as leptonic annihilation. Another sign could be in stud- 
ies of polarization. As observed in Ref. @ in the context 
of LO NRQCD, associated production leads to a more 
complex polarization structure than that expected from 
a pure octet mechanism. Although color transfer does 
not fall cleanly into the analysis of NRQCD, we might 
expect much the same spin symmetry to apply, so that 
a study of polarization effects can also help characterize 
this mechanism. 

To develop an effective theory that incorporates color 
transfer, we shall have to add new operators. As al- 
ready found in Ref. [9( in the context of NRQCD factor- 
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ization, the long-distance transition from octet to color 
singlet states can be associated with a nonlocal opera- 
tor of the form d\' A' [P^q" F v ^ a {\' P) ], with P and 
q as in Eq. ([1]). Noting that X'q can be interpreted as 
the configuration-space separation of a pair, we see that 
this is the same operator that represents the singlet-octet 
transition in the potential NRQCD Hamiltonian [HI, [H[, 
here extended to production. It is clear, however, that 
this correspondence holds only in the region where we 
can strongly order the velocities. More generally, we shall 
have to introduce products of three fields, as V^Xi^jv 
whose interactions with the soft gluon field will be labeled 
by momenta (as in NRQCD) according to heavy quark 
effective theory. In general the nonperturbative dynam- 
ics these operators will overlap with those of NRQCD, as 
well as those of potential NRQCD. It will require some 
care to avoid double counting and to develop a formalism 
that can be matched to the phenomenology of associated 
production. 
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